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DEFACE 


In  this  report  we  describe  the  individual  components  and  sequence  of 
instructions  of  a  computer  program  designed  for  the  st.ereological 
msasurement  of  fibers.  Utilizing  on-screen  menu  oanmards,  we  were  able  to 
design  and  execute  a  step-by-step  analysis  of  a  fiber  measurement 
program.  There  was  seme  difficulty  in  writing  the  program  due  to  a  lack 
of  clear  instructions  in  the  operator's  manual;  however,  conversations 
with  company  representatives  helped  us  overcome  these  problems. 


A  COMHJTER  EROGRAM  FOR  MEASURING  FIBERS 


WITH  THE  ZEISS  CSM  950  SCANNING  ELECTRON  MICROSOOPE 


IMIRODUCTTON 


The  abjective  of  this  report  is  to  describe  a  computer  program 
designed  for  measuring  stereological  features  of  fiber's  by  means  of  a 
Kontron  image  analysis  system,  which  is  built  into  a  Zeiss  CSM  950 
Scanning  Electron  Microscope  (SEM) . 

J  The  program  was  written  to  permit  interaction  between  the  menu  program 
displayed  on  the  IV  monitor  aivi  the  keyboard  or  mouse.  This  increased 
capability  permits  the  operator  to  change  parameters  within  the 
measurement  program  depending  on  the  type  of  fiber  sample. 

Since  research  presently  being  conducted  in  these  laboratories  on 
ballistic  fibers,  including  biotechnological ly "derived  fibers,  requires 

information  about  fiber  size,  shape,  etc.,  developing  a  menu*  lied 

!  I  . 

computer  program  facilitates  that  purpose.  J,  •  t  t  ■/*/  j .  -•  ^ 

MEIHODS 

The  sample  model  use  was  a  drawing  of  11  fibers  (Fig.  1)  which  were 
visualized  with  a  Dage  MTI  Model  CCD  725  video  camera  with  a  Fuji  C6  X 
17.5  TV  zoom  lens.  The  images  from  the  TV  camera  were  transferred  to  a 
Zeiss  CSM  950  SEM  with  a  built-in  Kontron  image  analysis  system.  The  CSM 
950  utilizes  a  CP/M  operating  system  having  a  Z80  processor  with  a 


FIGURE  1.  Drawing  of  II  Fibers  to  be  Measured. 

The  parameters  measured  (in  pixels) 
were  Area,  Perimeter,  Feret  X 
(maximum  projected  length  in  the 
horizontal  direction) ,  ard  DCircle 
(area  of  fibers  within  a  circle 
measuring  100  pixels  in  diameter) . 


2 
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20  Mbyte  hard  disk,  64  Kbyte  randan  access  memory  (RAM)  and  a  16  Kbyte 
video  RAM. 

The  inage  processing  system  performs  evaluations  in  a  series  of 
steps.  First,  the  image  fron  a  TV  camera,  SEM,  light  microscope  or  any 
other  internal  or  external  source  is  digitized  and  stored.  Next,  the  grey 
image  is  processed  for  improved  contrast  and  signal/noise  ratio.  Then, 
certain  features  are  extracted,  i.e. ,  unwanted  or  overlapping  features  are 
eliminated,  ard  features  such  as  size,  shape,  etc.,  are  analyzed. 

Finally,  the  derived  data  are  processed  statistically  and  a  printout  or 
photograph  is  obtained. 

There  are  several  major  functions,  which  are  .rouped  depending  upon 
the  specific  sequence  of  a  standard  measuring  price  dura.  The  functions 
axe; 

1.  Input  -  defines  the  type  of  input  device  and  controls  image 
digitization  and  storage. 

2.  Calibrate  -  controls  the  scaling  factors,  i.e.,  inches,  microns, 
etc. 

3.  Enhance  -  rescales  the  grey  levels,  thereby  enhancing  image 
quality.  Image  noise  reduction  fir'ters  are  used. 

4.  Image  Edit  -  contains  functions  for  the  editing  of  images. 

5.  Segment  -  extracts  the  image  background  (distinguishes  between 
objects  and  background) . 

6.  Multiphase  -  enables  further  processing  of  binary  (two  grey 
levels)  and  multiphase  (more  than  two  grey  levels)  images. 
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7.  Parameters  -  permits  the  selection  of  object-specific  (single 
objects  in  a  field  of  view) ,  or  field-specific  (area  of  a  visual  field 
minus  the  area  occupied  by  objects) . 

8.  Evaluate  -  evaluates  objects  and  outputs  the  results. 

9.  Advanced  -  performs  geometric  acid  arithmetic  transformation  of 
iirages,  i.e. ,  improves  the  quality  of  an  image  using  Fast  Fourier 
Transform  methods. 

10.  Utilities  -  has  general  auxiliary  functions. 

11.  Peripheral  -  permits  control  of  peripheral  devices,  remote 
control,  etc. 

All  the  above  functions  have  subfunctions  that  cure  displayed  on  the 
video  monitor  in  the  form  of  a  menu.  Dialogue  with  the  program,  allowing 
the  operator  to  assign  or  delete  numerical  values  as  well  as  to  control 
all  of  the  system  functions,  is  accomplished  by  utilizing  the  keyboard  or 
a  mouse. 

The  following  is  the  selection  sequence  for  the  operator  in  order  to 
obtain  the  desired  program. 

A.  Calibrate 

1,  Scale  -  As  all  measurement  data  are  indicated  in  the  selected  unit  of 
calibration,  the  most  inportant  prerequisite  for  every  geometric 
neasurement  is  calibrating  the  measuring  system.  Ihe  calibration 
procedure  is  carried  out  according  to  the  instructions  in  the  operator's 
manual  (Anonymous,  1986) . 

a.  SChQ,  the  number  of  the  scale  factor  to  be  activated,  is  set  at  1. 
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1.  Reset  -  Before  starting  a  new  measuring  sequence,  e.g. ,  after'  loading 
a  new  program,  sane  conditions  such  as  data  buffering,  identification 
number  and  scanning  stage  can  be  reset.  The  system  then  has  the  same 
status  as  after  initialization. 

a.  Single  (on)  -  Clears  the  buffer  into  which  all  the  data  from 
several  fields  are  accumulated. 

b.  Class  (on)  -  Clears  the  data  summed  up  over  several  evaluated 
fields  to  provide  average  values  (e.g. ,  mean  area  percentage) . 

C .  Select  Parameter 

1.  Area  -  Object  specific  area. 

a.  NdS  (20)  -  Number  of  classes  when  classification  is  selected. 
Classification  refers  to  a  range  of  measurement  parameter  values,  i.e. 
area,  perimeter,  etc.  into  which  classes  (groups)  of  fibers  waald  fit. 

b.  MODX  (1-3)  -  Presentation  and  scaling  mode  of  the  X-axis  in  the 
histogram  of  the  results. 

c.  MODY  (1-10)  -  Presentation  and  scaling  mode  of  the  Y-axis  in  the 
histogram  of  the  results. 

d.  Single  (on)  -  selection  of  a  single  list. 

e.  Class  (on)  -  Selection  of  a  classification. 

f.  Low  0.000  -  Lower  bound  of  classification. 

g.  High  5000  -  Upper  bound  of  classification. 

D.  Select  Parameter 

1.  Feret  -  Feret  diameters  (maximum  projected  lengths)  in  X  and  Y 
direction. 

a.  NCLS  (2.0)  -  Number  of  classes  when  classification  is  selected. 
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b.  MODX  (1-3)  -  Presentation  and  scaling  mode  of  the  X-axis  in  the 
histogram  of  the  results. 

c.  MODY  (1-10)  -  Presentation  and  scaling  mode  of  the  Y-axis  in  the 
histogram  of  the  results. 

d.  Single  (on)  -  Selection  of  a  single  list. 

e.  Class  (on)  -  Selection  of  a  classification. 

f.  Low  (0.000)  -  Lcwer  bound  of  classification. 

g.  High  (200.00)  -  Higher  bound  of  classification. 

h.  NCLY  (20)  -  Number  of  classes  when  classification  is  selected. 

i.  MDYX  (1-3)  -  Presentation  and  scaling  mode  of  the  X-axis  in  the 
histogram  of  the  results. 

j.  MDYY  (1-10)  -  Presentation  and  scaling  mode  of  the  Y-axis  in  the 
histogram  of  the  results. 

k.  LO.Y  (0.000)  -  Lower  bound  of  the  classification. 

l.  HI .Y  (1000.0)  -  Upper  bound  of  the  classification. 

E .  Select  Parameter 

1.  Perim  -  Object  specific  perimeter  and  length. 

a.  NCLS  (20)  -  Number  of  classes  when  classification  is  selected, 
b-  MODX  (1-3)  -  Presentation  and  scaling  mode  of  the  X-axis  in  the 
histogram  of  the  results. 

c.  MODY  (1-10)  -  Presentation  and  scaling  mode  of  the  Y-axis  in  the 
histogram,  of  the  results. 

d.  Single  (on)  -  Selection  of  a  single  list. 

e.  Class  (on)  -  Selection  of  a  classification. 
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f.  Low  (0.000)  -  Irwer  bound  of  classification. 

g.  High  (600.0)  -  Upper  bound  of  classification. 

h.  NCUT  (10)  -  Number  of  classes  when  classification  is  selected. 

i.  MDLX  (1-3)  -  Presentation  and  scaling  mode  of  the  X-axis  in  the 
histogram  of  the  results. 

j.  MDLY  (l-io)  -  Presentation  and  scaling  mode  of  the  Y-axis  in  the 
histogram  of  the  results. 

k.  LOUT  (600.0)  -  Lower  bound  of  classification. 

l.  HILT  (0.000)  -  Upper  bound  of  classification. 

F.  Select  Parameter 

1.  DCIRCL  -  Diameter  of  area-equivalent  circle  (area  of  fibers  within  a 
circle  which  measures  100  pixels  in  diameter) . 

a.  NCLS  (20)  -  Number  of  classes  when  classification  is  selected. 

b.  MODX  (1-3)  -  Presentation  and  scaling  mode  of  the  X-axis  in  the 
histogram  of  the  results. 

c.  MODY  (1-10)  -  Presentation  and  scaling  mode  of  the  Y-axis  in  the 
histogram  of  the  results. 

d.  Single  (on)  -  Selection  of  a  single  list. 

e.  class  (on)  -  Selection  of  a  classification. 

f.  Lew  (0.000)  -  lower  bound  of  classification. 

g.  High  (100.00)  -  Upper  bound  of  classification. 

G.  Select  Utilities 

1.  LAB:  Label  definition  -  Labels  are  used  as  markers  in  the  sequence  of 
a  measuring  program  to  which  a  jump  (analogous  to  a  GOTO  function)  can  be 
executed. 

a.  LEL# (1)  -  Identifier  of  the  label,  which  serves  as  designation  for 
a  jump. 
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1.  TVON  -  Switches  the  color  monitor  directly  to  a  TV  signal. 

a.  Online  (on)  -  Switches  the  online  status  on  or  off. 

I._  Select  Utility 

I.  SYNC  -  The  image  Processing  System  can  be  synchronized  externally 
(e„g. ,  by  using  a  TV  camera)  or  internally.  In  the  internal  mode,  no  TV 
input  is  possible. 

a.  Internal  (off) ,  Internal  (on) ,  or  External  (off) . 

J.  Select  Utility 

1.  Pause  -  Interrupts  run  of  program. 

K.  Input 

1.  TVINP  -  This  function  stores: 

a.  'iv  Input  -  A  TV  image  iji  the  memory  indicated  by  INT 

b.  INP  (1)  -  Image  memory  into  which  the  image  is  to  be  stored. 

L-  Choose  Segment. 

1.  DISC2L  -  This  function  separates  objects  from  the  background,  by 
setting  two  thresholds.  Either  the  grey  levels  inside  or  those  cutside 
the  entered  limits  are  set  to  grey  value  0  (black)  and  constitute  the 
background.  Depending  on  the  variable  BINARY ,  the  remaining  object  points 
either  keep  their  original  grey  levels  or  are  set  to  white  (255) . 

a.  INP  (1)  -  Image  memory  containing  the  image  to  be  processed. 

b.  OUT  (2)  -  Image  memory  into  which  the  discriminated  image  is 
stored. 

c.  LEVI.  (0)  -  lower  discrimination  limit, 
ri.  LEV2  (118)-  Upper  discrimination  limit. 
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e.  BINARY  (on)  -  The  object  pixels  become  255  (white) ,  i.e.  a  binary 
image  is  produced. 

M.  Calibrate 

1.  MFRAME  -  This  function  is  used  to  create  a  rectangular  or  circular 
measuring  frame,  which  is  sometimes  necessary  to  correct  errors  (eliminate 
overlapping  images)  caused  by  the  finite  size  of  the  field  of  view  during 
the  measurement  process. 

a.  INP  (2)  -  Memory  containing  the  image  to  be  processed. 

b.  WDSX  (508)  -  X-side  of  the  rectangular  frame. 

c.  WDSY  (467)  -  Y-side  of  the  rectangular  frame. 

d.  XO  (4)  -  X  position  of  the  (rectangular  or  circular)  frame. 

e.  YO  (6)  -  Y  position  of  the  (rectangular  or  circular)  frame. 

f.  RAD  (100)  -  Radius  of  the  circular  frame. 

N.  Evaluate 

1.  IDEOT  -  This  function  serves  two  purposes:  Identification  of 
discriminated  images  and  elimination  of  objects  following  the  conditions 
set  by  the  measuring  frame  (function  MFRAME,  group  Calibrate) . 

a.  INP  (2)  -  Image  memory  containing  the  image  to  be  identified. 

b.  OJT  (3)  -  Memory  containing  the  identified  and  frame-corrected 
image. 

c.  MARG  (1)  -  Controls  the  meaning  of  the  measuring  frame  defined  by 
MFRAME  (function  group  Calibrate) .  This  determines  whether  or  not  the 
fiber  touching  the  rectangular  or  circular  frame  on  the  video  monitor 
is  measured.  For  example,  one  might  not  want  to  measure  a  fi-ber 

ext ending  beyond  the  frame's  margin. 
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d.  8 -CONN  (on)  -  Method  of  identi  f  icaticn,  on  =  8  ODNN. 


...  Evaluate 

1.  MEASURE  -  This  function  executes.-  the  measurement  of  selected 
parameters  such  sis  area,  perimeter,  etc. 

a.  OBJ  -  Must  be  switched  on  if  object-specific  parameters  are  among 
the  previously  selected  parameters  (AREA-TIME) . 

b.  INP  (3)  -  Image  memory  containing  the  discriminated  OBJ  image. 

c.  GRIM  (1)  -  Image  memory  containing  the  original  grey  image. 

d.  AUXl  (9)  -  Auxiliary  image  memory. 

e.  ALJX2  (9)  -  Auxiliary  image  memory. 

f.  CHAN  (1)  -  Channel  selection. 

g.  SPEC  (1)  -  Distance  in  pixels  between  consecutive  lines  of  the 
"line  grid".  The  lines  are  simply  lines  of  the  image  memory,  or  pixel 
addresses  with  varying  X,  but  constant  v  coordinates. 

h.  FIG#  (1)  -  logical  flag  indicating  a  buffer  overflow. 

i.  REJECT  (on)  ~  Interactive  rejection  of  undesired  objects, 
preceding  an  automatic  measurement. 

P.  Evaluate 

1.  CUTCLS  -  (Output  of  Classification)  This  function  displays  the 
results  of  the  classification  in  the  form  of  a  list  and/or  as  a  histogram. 

a.  HISTO  (on)  -  Provides  measured  data  in  histogram  form. 

b.  HALT  (on)  -  Presents  data  on  the  monitor;  the  sum  histogram  can  be 
manipulated. 
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c.  NFAC  (0.000)  -  Defines  the  area  (in  user  units)  into  which  the 
absolute  counts  should  be  normalized.  A  value  of  0.0  gives  normalized 
absolute  frequencies. 

0.  Evaluate 

1.  OOTSGL  -  (Output  of  single  list)  This  function  displays  the 

individual  lists  of  all  measured  data. 

a.  HALT  (on)  -  Displays  lists  on  the  data  monitor  after  each 
measurement. 

R.  Select  Utility 

1.  Pause  -  This  function  interrupts  the  running  of  a  measuring  program. 
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RESUIIIS 


The  computer  program  designed  and  written  for  the  stereological  analysis 
of  several  fibers  (Table  1)  was  utilized  to  measure  their  Area,  Perimeter, 
Fenet  X  and  D  Circle.  The  inputs  to  the  parameters  of  this  program  were  as 
follows: 

Area  (high  upper  bound  of  classification)  -  5000 
Feret  X  (higher  bound  of  classification)  -  200 
Feret  X  (upper  bound  of  classification)  -  1000 
Perimeter  (high  upper  bound  of  classification)  -  600 
DCircle  (high  upper  bound  of  classification)  -  100 

The  printout  of  the  computer  program  can  be  seen  in  Table  1,  the  plots  for 
each  parameter  can  be  seen  in  Figures  2  to  5  and  the  printout  of  the 
statistical  data  can  be  seen  in  Table  2. 


CONCLUSIONS 

A  computer  prcgram  was  designed  and  written  to  measure  certain 
stereological  parameters  of  fibers  and  to  generate  subsequent  statistical 
analysis  of  the  results.  The  program  was  developed  to  be  driven  by  a  series 
of  menu  commands,  which  allowed  for  direct  interaction  between  the  image 
appearing  on  the  video  monitor  and  the  keyboard. 


TABLE  1.  Image  Analysis  program  for  Measuring  Fibers  Using  the  Zeiss  CSM 
950  SEM  with  a  Built-in  Kontron  Image  Analyzer. 


MAJOR 

FUNCTIONS 


SUB  INTEGER  VARIABLES  SCALE 

FUNCTIONS  VALUES  VALUES 


A. 

CALIBRATE 

SCALE 

SCNO 

1 

B. 

EVALUATE 

RESET 

SINGLE 

CLASS 

PARAM 

C. 

PARAMETERS 

AREA 

NCLS 

20 

SINGLE 

LOv 

0.000 

MODX 

1 

CLASS 

HIGH 

5000.0 

MODY 

1 

D. 

PARAMETERS 

FERET 

NCLS 

20 

SINGLE 

LOW 

0.000 

MODX 

1 

CLASS 

HIGH 

200.0 

MODY 

1 

LO.Y 

0.000 

NCLY 

20 

HI  -Y 

1000.0 

MDYX 

1 

MDYY 

1 

E. 

PARAMETERS 

FERIM 

NCLS 

20 

SINGLE 

LOW 

0.000 

MODX 

1 

CLASS 

HIGH 

600.0 

MODY 

3. 

LOUT 

600.0 

NCLT 

10 

LILT 

0.000 

MDLX 

1 

MDLY 

1 

F. 

PARAMETERS 

PCIRCL 

NCLS 

20 

SINGLE 

LOW 

0.000 

MODX 

1 

CLASS 

HIGH 

100-0 

.  MODY 

1 

G. 

UTILITIES 

LABI 

LBL* 

1 

H. 

iNRjr 

TVON 

ONLINE 

I. 

UTILITIES 

SYNC 
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(Continued) 


TABLE  1.  Image  Analysis  Program  for  Measuring  Fibers  Using  the  Zeiss  CSM 
950  SEM  with  a  Built-in  Kontron  Image  Analyzer.  (Continued) 


MAJOR 

SUB 

INTE3GER 

VARIABLES 

SCALE 

FUNCTIONS 

FUNCTIONS 

VALUES 

VALUES 

J.  UTILITIES 

PAUSE 

K.  INPUT 

TVINP 

INP  1 

L.  SEGMENT 

DISC2L 

LEVI  0 

IEV2  118 

BINARY 

LEVI  0 

M.  CALIBRATE 

MFRAME 

INP  2 

WQSX  508 

WDSY  467 

X0  4 

YO  6 

RAD  100 

N.  EVALUATE 

IDENT 

INP  2 

CUT  3 

MARG  1 

8-CONN 

0.  EVALUATE 

MEASUR 

INP  3 

GRIM  1 

AUX1  9 

OBJ 

AUX2  9 

CHAN  1 

SPAC  1 

FLGs 

REJECT 

P.  EVALUATE 

OUTCLS 

HIST 


NFAC  0.000 


Q.  EVALUATE  OOTSGL 


HALT 


R.  tmUTIES  PAUSE 
TOTAL  PROGRAM  LENGTH 


340  BYTES 


HALT 


14 


AES.  FREQUENCY 


COUNTS 

11 

UND.  FLOW 

0 

OVERFLOW 

0 

CLASSES 

20 

%-IL  CLASSES 

25% 

7 

MEDIAN 

9 

75% 

11 

95% 

18 

WDDUL 

250.0 
L.  BOUND 
0.000 
U.  BOUND 
5000. 
V  1~KT  1M 

1142. 

MAXIMUM 

4254. 

MEAN 

2177. 

SD 

901.3 

AREA 


FIGURE  2.  Histogram  of  Area  Measurements  from  Figure  1. 

The  X  axis  values  are  in  pixels. 

The  Y  axis  represents  the  frequency  disti  lbution . 
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ABS.  FREQUENCY 


10.0 


8.00 


6.00 


j  4.00 

2.00 


0.00 


COUNIS 

11 

UND.  FLOW 

0 

OVERFLOW 

0 

CLASSES 

20 

%-IL  CLASSES 

25% 

8 

MEDIAN 

9 

75% 

11 

95% 

19 

MDDUL 

30.00 
L.  BOUND 
0.000 
U.  BOUND 
600.0 
MINT!’ DM 
16V.  7 
rtfOCTHBN 

547.3 
MEAN 

289.4 
SD 

108.2 

PERIM 


0.00  200.0 


400.0  600.00 


FIGURE  3.  Histogram  of  Perimeter  Measurements  from  Figure  1. 
The  X  axis  values  are  in  pixels. 

The  Y  axis  represents  the  frequency  distribution. 
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ABS.  FRBQtOCY 

10.0  - 


COUNTS 
UND.  FLOW 
OVERFLOW 
CLASSES 
*-IL  CLASSES 
25% 

MEDIAN 


8.00  -4 


75% 

95% 
MODUL 
20.000 
L.  BOUND 
0.000 
U.  BOUND 
200.0 
MINIMUM 

■*  r\r\ 
Jli  vw 

MAXIMUM 

191.0 

MEAN 

77.91 

SD 

48.38 
FERET  X 


11 

0 

0 

20 


FIGURE  4.  Histogram  of  Feret  X  Measurements  firm  Figure  1. 
Ihe  X  axis  values  are  in  pixels. 

Hie  Y  axis  represents  the  frequency  distribution. 
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m  \o  «-i  o 

ri  CN 


ABS. 

10.0 

8.00 

6.00 

4.00 

2.00 

0.00 


FREQUENCY 


*  ■  ■  ■  1  t  1  1 


COUNTS 

11 

UND.  FLOW 

0 

OVERflCW 

0 

CLASSES 

20 

%-IL  CLASSES 

25% 

9 

MEDIAN 

11 

75% 

12 

95% 

15 

MDDUL 

L. 


5.000 
BOUND 
0.000 
U.  BOUND 
100.00 
MINIMUM 

-»C3  11 

MAXIMUM 

73.60 

MEAN 

51.71 

SD 

30.37 

-DCIRCLE 


0.000 


20.00 


40.00 


60.00 


80.00 


100.0 


FIGURE  5.  Histogram  of  an  Area  Equivalent  to  a  Circle 
(called  a  DCircle)  within  Figure  1. 

The  X  axis  values  are  in  pixels. 

The  Y  axis  represents  the  frequency  distribution. 
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TABLE  2.  Statistical  Data  (in  Pixels)  for  11  Fibers  Measured 


COUNT 

AREA 

PERIM. 

FERETX 

DCTRCLE 

1 

3111.0 

413.4 

137.0 

62.94 

2 

2134.0 

270.0 

53.0 

52.13 

3 

2512.0 

318.4 

58.0 

56.55 

4 

1731.0 

250.8 

103.0 

46.95 

5 

2378.0 

302.6 

52.0 

55.03 

6 

1691.0 

226.0 

41.0 

46.40 

7 

2196.0 

269.5 

44.0 

52.88 

8 

1142.0 

169.7 

68.0 

38.13 

9 

1217.0 

196.7 

79.0 

39.36 

10 

1582.0 

219.2 

31.0 

44.88 

11 

4254.0 

547.3 

191.0 

73.60 
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